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Abstract 

Many habitat specialist species are originally composed of small, discontinuous populations because their habitats are 
naturally fragmented or patchy. They may have suffered the long-term effects of natural patchiness. Mediterranean 
heathlands, a representative habitat in the Strait of Gibraltar region, are associated with nutrient-poor, acidic sandstone 
soils. Sandstone soil patches in the African side of the Strait (Tangier) are, in general, smaller and more scattered than in the 
European side (Algeciras). In this study, we analyze the effect of this sandstone patchiness on the population genetic 
diversity and structure of two Erica species from these Mediterranean heathlands that differ in their edaphic specificity, £ 
australis, sandstone specialist, and £ arborea, generalist. Average levels of within-population genetic diversity and gene flow 
between populations were significantly lower in Tangier (high sandstone patchiness) than in Algeciras (low patchiness) for 
the sandstone specialist, whereas no differences between both sides of the Strait were detected in the edaphic generalist. 
Since most endemic species in Mediterranean heathlands of the Strait of Gibraltar are sandstone specialists, these results 
highlight an increased vulnerability to loss of genetic diversity and local extinction of the heathland endemic flora in the 
Tangier side of the Strait of Gibraltar. 
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Introduction 

Plant populations in small habitat fragments have frequently 
smaller population sizes and experience a higher degree of 
isolation among populations than those from continuous or less 
fragmented habitats [1]. The smaller and more isolated the 
populations, the more liable they are to demographic stochasticity, 
to genetic depletion and, ultimately, to local extinction [2—4]. 
These small, fragmented populations are also highly vulnerable to 
edge effects, habitat degradation or to catastrophic disturbance 
events, either natural or human-induced [5,6]. 

Genetic diversity in a plant population is lost mostly by the effect 
of genetic drift and by reduced genetic communication (i.e. low 
gene flow and migration) among populations [7]. The loss of 
genetic diversity negatively affects the adaptability of populations 
to environmental changes as it constitutes the basic evolutionary 
substrate and therefore lessens their chances of survival [8,9]. 
Gene flow reduces the risk of extinction of small populations either 
by avoiding genetic erosion through the input of new alleles mainly 
via seed dispersal or by recolonization after a local extinction event 
[10— 1 1]. Therefore, gene flow allows the long-term persistence of 
small populations and it is a critical process for the management 
and conservation of threatened species in fragmented or degraded 
habitats [3,12]. 

Habitat degradation and fragmentation have indeed become 
major research topics in both population genetics and conserva- 



tion biology [13], and many scientific contributions report the 
effects of anthropogenic fragmentation on the population genetic 
structure of species that previously occupied large, continuous 
habitats (e.g. [3,14,15]). However, fragmentation has not always a 
human-induced cause: many habitat specialist species are origi- 
nally composed of small, discontinuous populations because their 
habitats are naturally fragmented or patchy [16]. Those habitat 
specialist species may have undergone the long-term effects of 
natural patchiness and their populations might exhibit low genetic 
diversity within populations and strong genetic structure among 
populations (e.g. [17]; but see [18]). As a result, a specialist species 
confined to a naturally patchy habitat would be inherendy more 
sensitive to further fragmentation and/ or degradation of the 
habitat than a generalist species occurring within and outside those 
habitat patches [19,20]. It is thus crucial to ascertain whether 
natural habitat patchiness does actually affect the genetic diversity 
of populations of habitat specialist species. 

Many endemic, and/ or endangered plant species in the 
Mediterranean Basin, one of the world's biodiversity hotspots, 
are habitat specialists — edaphic and/or orographic — restricted 
to uncommon, patchy habitats. Mediterranean heathlands of the 
Strait of Gibraltar region, at the western end of the Mediterra- 
nean, harbour an endemic-rich flora composed mainly of edaphic 
specialist species, restricted to nutrient-poor, highly acidic, 
sandstone soils (sandstone specialists; [21,22]). Mediterranean 
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heathlands on either side of the Strait of Gibraltar are very similar 
from both floristic and landscape points of view [23]. However, the 
absolute number and relative abundance of sandstone endemics in 
heathland communities of the African side of the Strait (hereafter 
Tangier) are lower than those from the European side (hereafter 
Algeciras) [22]. This seems to be a direct outcome of a striking 
feature of the Strait of Gibraltar: the generally smaller size and 
more isolated pattern of the sandstone patches in Tangier (Fig. 1 ; 
[22]). The overall higher sandstone patchiness in Tangier might 
limit gene flow and recolonization among populations of sandstone 
specialists — most of them endemics — thus making them highly 
prone to local extinction [19] in comparison with the generally 
more continuous and larger sandstone patch sizes in Algeciras 
(Fig. 1). 

Here we explore whether an overall high level of sandstone 
patchiness is associated to low levels of within-population genetic 
diversity and of gene flow between populations of sandstone 
specialist species. To do so, we studied the genetic diversity and 
population structure of two heath species, Erica australis L. and E. 
arborea L. by means of microsatellite markers. These two species 
are abundant in Mediterranean heathlands from both Algeciras 
and Tangier, but they differ notably in their edaphic (sandstone) 
specificity: while E. australis is a sandstone specialist, virtually 
restricted to sandstone-derived soils on both sides of the Strait, E. 
arborea has a broader edaphic spectrum [24]. We compared 
population genetic diversity indices, including relative levels of 
gene flow among populations, for each species between the two 
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Figure 1 . Geographical location of sampling sites of £ australis 
and £ arborea. At each site in Algeciras (black dots) and Tangier (white 
dots), one population of each species was sampled (except in A14, grey 
dot, where only E. australis was sampled, see Table 1). Grey shaded 
patches indicate the presence of nutrient-poor, highly acidic, sandstone 
soils. 

doi:1 0.1 371 /journal.pone.0098602.g001 



sides of the Strait as surrogates for overall low (Algeciras) and high 
(Tangier) levels of sandstone patchiness. Since an overall high 
patchiness would affect sandstone specialists more strongly than 
generalists, we expect to find lower levels of genetic diversity and 
gene flow in Tangier than in Algeciras populations for E. australis, 
but not for E. arborea. 

This study provides valuable insights into the differential effects 
of natural habitat (sandstone) patchiness on the genetic diversity 
and structure of populations of habitat specialist and generalist 
species. Considering that most endemic species from these 
Mediterranean heathlands are sandstone specialists [21,22,25], 
this study also allows us to highlight a presumably higher 
vulnerability to loss of genetic diversity and ultimately local 
extinction of the heathland endemic flora in the African side of the 
Strait of Gibraltar (i.e. Tangier), probably as a consequence of its 
overall higher sandstone patchiness. From a conservation 
perspective, this knowledge could help focus appropriate conser- 
vation on those patchy communities to better maintain the high 
biodiversity and uniqueness of Mediterranean heathlands across 
the Strait of Gibraltar region. 

Materials and Methods 

Ethic Statement 

Necessary permits for fieldwork and sampling were obtained 
through the University Abdelmalek-Essaadi (Morocco) and the 
Junta deAndalucia (Spain). None of the two species used in this study 
is endangered or red-listed. 

Study Species 

Erica australis and E. arborea (Ericaceae) are two relatively 
abundant heath species in Mediterranean heathlands on acid, 
nutrient-poor, Oligo-Miocene sandstone soils at both sides of the 
Strait of Gibraltar [25]. They are both diploid (2«= 2x = 24; [26]) 
and monophyletic [27,28], and may be considered genetically 
related as inferred from a high transferability rate from SRR 
markers from a South African congener [29]. They are also 
morphologically and ecologically similar (e.g. long-lived, woody 
plants with ericoid leaves, acidophilous and post-fire resprouters; 
[24,30,31]. However, they differ notably in their edaphic 
specificity: E. australis is a sandstone specialist, restricted to highly 
acidic, sandstone soils, whereas E. arborea has a broader edaphic 
spectrum. Within the regional context of this study (i.e. Strait of 
Gibraltar region), E. australis is virtually restricted to highly acidic, 
aluminium-rich sandstone soil patches whereas E. arborea occurs 
both within sandstone patches and out in the surrounding matrix 
on different, non-acid soils [24] . They also differ in the amplitude 
of their geographical range: E. australis is endemic to the western 
end of the Mediterranean basin (western third of the Iberian 
Peninsula and NW tip of Africa), whereas E. arborea has a much 
wider, circum-Mediterranean and eastern African distribution 
[27,30]. 

Regarding aspects of reproductive biology, both species have 
small, hermaphroditic flowers, pink-coloured in E. australis and 
whitish in E. arborea. They are insect pollinated, with pollen shed in 
tetrads [32] . Erica australis sets flowers in late-winter to early-spring 
(December to April), whereas E. arborea is a spring-blooming 
species (March to May; [33]). Both species produce small, oval 
seeds that lack any kind of appendages for long distance dispersal. 
Erica australis seeds are larger (0.9-1.1 mm length) than those of E. 
arborea (0.4-0.5 mm; [34]). Such small seeds may be wind- 
dispersed but only over limited distances (less than 100 m; [35]). 
Erica australis presents a small caruncle [32], which might indicate 
also short-distance, ant-dispersal. 



PLOS ONE | www.plosone.org 



2 



May 2014 | Volume 9 | Issue 5 | e98602 



Patchiness and Population Genetic Structure 



! i 



00 \£> CTi Oi CO CO 



O t- «- 



3 5 



S S ^ 



vo m vo 



mmmLnLnvOLnvo 



•— \o in 



o i— in i— 



(N lo m vo 



mu-iuni-nvOLnLnm 



\o vo vo 



m vo m m vo vo vo 



ro ro «— 



ro ro ro ro ro ro 



I ! ! i ! ! 



ill! 



fN Oi 



CO vO m Ov 



! ! ! 



d m ^ t rN \D 

i- i- o 



«- o 

CO 

O r~ 



r 1/1 CO 



1° 



CT\CXi< — in ' — LnoorNCiuim 
lOinvOvOvOvDvOvOvOvOvOvO 



\0 vO vO 



' — >— CO C* CO 



ro (N ro ro 



r^^-corMLorSi— ^-o 
d> ^ d d d> 

t— CTi CTi CO r- vO >— CO*— 



mrorororNrNroro 



555555555555 



5 5 5 



5 5 5 5 5 5 5 



ov i- i- 



ui lo m m 



IN t >- O (N 

Lfl LTl (N IN (N 



m in m m 



vo vo vo *o 



fN vO 



ro ro ro ro 



mmmmmmmm 
rorororororororo 



. o 

aJ ^ 

I s 

£ vO 



2 O Q) 



a 
o 



o 

< < 



www 
* 10 

o o 

< < 



13 

0 
+j 

<5 

3 

a. 

oo 
o 
< 



0 ^ 

OA A 

O i- 

< <C 



< 



o 
E 

3 

2 x 

01 TJ 

2 & 

a to 
-i a. 

m * 

3 < 



E 
o 



o- c 

■g * 3 £ 

" tj o 

Q C ED X 

£ * — — 

a h. ll oj m 

■- ai — .q .o 

oi "¥ 7 7 

N c i- in >o 

*- m o o o 

< I- I- I- H 



E 

01 

I 1 

CT1 0) 

* j= 

m vo 



a ^> 
£ 2 



PLOS ONE I www.plosone.org 



3 



May 2014 | Volume 9 | Issue 5 | e98602 



Patchiness and Population Genetic Structure 



Population Sampling, DNA Extraction and Microsatellite 
Amplification 

Twenty-two localities (11 from Algeciras and 11 from Tangier) 
where E. australis and E. arborea co-occurred, plus one locality in 
Algeciras where only E. australis was found, were sampled (Table 1 , 
Fig. 1). Population sampling consisted of fresh leaves from up to 30 
individuals per population (spaced at least 10 m from each other), 
except in four populations of E. australis and three populations of 
E. arborea where this number could not be reached (Table 1). Thus, 
the total sampling included 676 individuals of E. australis (357 from 
Algeciras and 319 from Tangier) and 631 individuals of E. arborea 
(308 from Algeciras and 323 from Tangier). Samples were dried in 
silica gel and stored at room temperature until DNA extraction. 
Dry leaf material, approximately 1 00 mg per sample, was reduced 
to fine powder using stainless steel beads on a Mixer Mill MM400 
cell disrupter (Retsch, Llanera, Spain). DNA was extracted using 
SpeedTools plant DNA extraction kit (Biotools, Madrid, Spain), 
and eluted in 50 (Xl in Tris-EDTA 0.1 xbuffer. 

Amplification of microsatellite loci followed Segarra-Moragues 
et al. [29,36]. Seven microsatellite loci (Ecocl08, Ecocll7, 
Ecocl32, Ecocl37, Ecocl42, Ecoc431 and Ecoc446) were 
amplified in E. australis, and eight (the same as for E. australis 
plus Ecoc 1 1 5) in E. arborea. These microsatellite loci were unlinked 
[29] and polymorphic in both species. Fifteen percent of the 1307 
genotyped individuals were included as duplicates to check for 
possible genotyping errors. All of them showed identical allelic 
profiles to their corresponding original samples. 

PCR products were electrophoresed in an ABI3730 automated 
sequencer (Applied Biosystems, Madrid, Spain) using LIZ500 as 
internal lane size standard. Assignment of fragments to allele 
classes was carried out with Genemarker version 1.97 software 
(Softgenetics, State College, PA). Genotypic matrices were 
deposited at DRYAD (http:/ /datadryad.org/) under accession 
number doi:l 0.506 1/dryad.bj 70k. 

Population Genetic Analyses 

GENETIX v. 4.05 [37] was used to estimate allele frequencies, 
mean number of alleles per locus (A), and observed (H 0 ) and 
unbiased expected (H E ) heterozygosities [38]. Wright's F-statistics 
were estimated according to Weir & Cockerham [39] using 
GENEPOP'007 [40] and tested for significance using Fisher's 
exact tests. This latter software was also used to estimate the 
Maximum likelihood frequency of null alleles and 95% confidence 
intervals for null allele frequencies according to Dempster et al. 
[41]. 

To evaluate, for each species, whether population diversity 
indices differed between Algeciras and Tangier groups of 
populations, average allelic richness per locus (A ) estimated 
according to the rarefaction method of Hurlbert [42] adapted by 
El Mousadik & Petit [43], average observed heterozygosity {H 0 ), 
average genetic diversity within populations (H s ), and inbreeding 
coefficient (.Fis), were compared using FSTAT version 2.9.3.2 [44] 
and tested for significance using 10,000 permutations. This same 
software was used to check for significant differences in population 
differentiation (i.e. average F ST values) for each species between 
Algeciras and Tangier. 

Bayesian analyses in STRUCTURE v. 2.1 [45,46] were used to 
estimate the population genetic structure and to infer the most 
likely number of genetic clusters (K) in both heath species. 
STRUCTURE assigns individuals to the K different genetic 
clusters based on allele frequencies at each locus. Estimated 
number of A" clusters ranged from 1 to 21, and analyses were based 
on an admixture ancestral model with correlated allele frequen- 



cies. In each run of the program Monte Carlo Markov Chain 
(MCMC) and burn-in period length consisted of 1.2 xlO 6 and 
8 x 1 0 5 iterations, respectively. The amount of variation of the 
likelihood was evaluated by carrying out ten runs for each K. The 
most likely number of genetic clusters (K) was estimated following 
Evanno et al. [47], which uses an ad hoc parameter (AK) to estimate 
the rate of change of likelihood values between successive iT values. 
POPULATIONS version 1.2.3. beta [48] was used to compute 
pairwise D A genetic distances [49] between populations of each 
species. These D A genetic distance matrices were used to compute 
eigenvalues and eigenvectors to perform Principal Coordinates 
Analyses (PCoA) for each species. Minimum Spanning Trees 
(MST) using D A distance matrices were constructed with 
NTSYSpc version 2.1 [50] and were superimposed onto the 2D- 
PCoA plots. MST computes minimum-length pairwise connec- 
tions between points. When superimposed onto PCoA plots it 
helps detecting local distortions (i.e. pairs of points which look 
close together in the PCoA but actually are far apart if other 
dimensions). Isolation by distance (IBD) was estimated separately 
for E. arborea and E. australis at each side of the Strait by matrix 
correlation analyses using a matrix of log-transformed pairwise 
geographical distances between populations and a matrix of a 
pairwise linearized F S t values (i.e., F ST /(l - F & t); [51]) computed 
with ARLEQUIN v. 3.1 1 [52]. Significance of the correlation was 
tested for each species and side of the Strait with Mantel test (1000 
permutations) using NTSYSpc. Then, ANCOVAs were used to 
test for parallelism of the resulting IBD slopes in Algeciras and 
Tangier within each species. Finally, gene flow between popula- 
tion pairs was also estimated separately for each species at each 
side of the Strait as the effective number of migrants per 
generation (Mm), obtained directly from F ST values by using 
Wright's [53] island model. Although these values should not be 
taken as accurate estimates of numbers of migrants, the 
comparison of average jV m values between populations at each 
of the two sides of the Strait may still be useful to explore for each 
species the existence of a decrease in gene flow associated with 
higher patchiness. For the sake of this study, we focused on 
geographical distances separating pairs of populations smaller than 
30 km in order to factor out or minimize possible geographical 
effects other than sandstone patchiness. 

Results 

Genetic Variation in E. australis and E. arborea 
Populations 

The seven microsatellite loci amplified a total of 167 alleles in 
Erica australis (Table SI). Number of alleles per locus ranged from 
10 (Ecoc446) to 41 (Ecoc 117), with a mean of 23.85± 10.02 alleles 
per locus. Mean number of alleles per locus ranged from a 
minimum of 6.42 from the Tangier population T12 to a maximum 
of 11.14 in the Algeciras populations A02 and A15 (Table SI). 
Average observed heterozygosities (Ho) ranged from 0.505 
(population T16) to 0.704 (population T06), and unbiased 
expected heterozygosities (H E ) ranged from 0.554 (population 
T16) to 0.786 (population A02; Table 1). In E. arborea a total of 108 
alleles were scored from eight microsatellite loci (Table S2). 
Number of alleles per locus ranged from 2 (Ecoc 11 7) to 30 
(Ecoc446) with an average of 1 3.50±9. 16 alleles. Mean number of 
alleles per population in E. arborea ranged from a minimum of 5.37 
in the Tangier populations T06 and T18 to a maximum of 7.37 in 
the Tangier population T16. Average observed heterozygosities 
ranged from 0.520 (population A09) to 0.650 (population A17), 
and unbiased expected heterozygosities ranged from 0.534(popu- 
lation T06) to 0.623 (populations A08 and A12; Table 1). 
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Significant deviation towards heterozygote deficiency (i.e. 
positive Fjs) was found in 22 out of the 23 populations studied 
of the sandstone specialist E. australis (Table 1). A similar pattern 
was also detected in the edaphic generalist E. arborea, although 
positive FIS values in this species were significant only in 1 0 out of 
the 22 populations (Table 1). High ,F IS values may be influenced by 
biological and methodological factors including inbreeding, local 
population substructure (i.e wahlund effect) and null alleles. 
Estimation of null allele frequencies gave significant results (i.e. 
lower 95%CI allele frequency >0.05) in 24 out of 161 and 39 out 
of 176 locus x population combinations in E. australis and E. 
arborea, respectively. Most of these significant combinations were 
concentrated in one locus (Ecoc431, 15 populations) in E. australis 
and two loci (Ecocl32, 22 populations and Ecoc431, 15 
populations) in E. arborea. Although the actual presence of null 
alleles in the dataset may only be confirmed by progeny analyses, 
our large sampling sizes provided indirect evidence that the 
frequencies of null alleles could have been overestimated by 
Dempster et al.'s (1977) method used in GENEPOP. As an 
example, in locus Ecoc431 we should have found ca. 20 null 
homozygotes out of 446 individuals genotyped in 1 5 populations 
that showed significant estimated frequencies of null alleles, but 
none was found. Therefore, null alleles, if at all present, should 
have litde impact in our datatset. Besides, repeated estimation of 
genetic diversity indices and population relationships did not 
change significandy after removing the two loci (data not shown). 
Accordingly, high Fi$ in these heath species may be related to 
pollen dispersal system favouring inbreeding and/or to population 
substructure, as it has been reported in other Erica species [36,54]. 

Although average within-population genetic diversity indices 
(A*, H 0 and H s ) in both species were slighdy lower in Tangier, 
where sandstone patchiness is more prominent (Fig. 1), only 
average H s was significantly lower in Tangier for the sandstone 
specialist E. australis, (Table 2). In this species, Algeciras 
populations showed a significantly higher inbreeding coefficient 
(Fjs) than Tangier ones, whereas no significant differences were 
detected in the edaphic generalist E. arborea (Table 2). 

Population Genetic Structure off. arborea and £ australis 

Pairwise F ST values in E. australis ranged from 0.010 (pair A02- 
A15) to 0.164 (T05-T16), whereas for E. arborea they ranged from 
zero (pair T01-T05) to 0.106 (pair T06-T16). Average F ST values 
were significandy different from zero in the two species, but higher 
in the sandstone specialist E. australis (average -Fst = 0.074) than in 
the edaphic generalist E. arborea (average F ST = 0.039). The 



comparison of average F ST values between the two sides of the 
Strait revealed for the two species a slightly stronger genetic 
differentiation in Tangier populations, although these differences 
were not significant in either species (Table 2). 

The method of Evanno et al. [47] for the estimation of the most 
likely number of genetic clusters (K) revealed a maximum modal 
value of AATor K—1 for both E. australis and E. arborea (see Figure 
SI). These values were AA"= 1794.02 and AK= 24.43 for E. 
australis and E. arborea, respectively (Fig. SI). A second maximum 
modal value of AA'=263.17 was found for A'=4 for E. australis 
(Fig. SI). The proportion of membership of the populations to 
each of the two genetic clusters (K— 2) in E. australis did not 
correspond to the geographical membership of the populations, 
with some populations of each side showing a higher proportion of 
membership to cluster 1 and others to cluster 2 (Fig. 2a). For K— 4, 
populations of E. australis from Algeciras showed a higher 
proportion of membership to clusters 3 and 4, whereas populations 
from Tangier showed a higher proportion of membership to 
clusters 1, 2 and 4 (Fig. 2b). Interestingly, populations of the 
sandstone specialist E. australis from the high sandstone patchiness 
area (Tangier) showed higher proportion of membership to more 
different genetic clusters than populations from the low patchiness 
one (Algeciras). On the other hand, in the generalist E. arborea, 
populations and individuals showed an almost symmetrical 
membership to both genetic clusters revealed at K= 2 in either 
Tangier or Algeciras (Fig. 2c). The method of Evanno et al. [47] is 
computationally constrained to detect K>2 genetic clusters; the 
most likely outcome for E. arborea is in fact K— 1 [55] which cannot 
be tested with this method. However, this homogeneous admixture 
pattern altogether with the low AK value obtained in E. arborea is 
indicative of an absence of population genetic structure (i.e. K= 1) 
throughout the Strait of Gibraltar and within Tangier and 
Algeciras [45] suggesting a weaker effect of patchiness in this 
species (Fig. 2c). 

The two-dimensional Principal Coordinates Analyses (PCoA)+ 
MST revealed similar results to those from STRUCTURE, and 
also showed direct links between populations (Figs. 3a-b). In E. 
australis the PCoA separated two population groups along the first 
axis, both of them including Algeciras and Tangier populations 
(Fig. 3a). In E. arborea, by contrast, the spatial arrangement of the 
populations in the PCoA was somewhat consistent with their 
geographical membership (Fig. 3b). The largest group was 
composed of two subclusters of populations. One of these included 
all Algeciras populations, with the exception of population A09, 
which appeared derived from the second subcluster that included 
all northernmost coastal Tangier populations. This second 



Table 2. Comparison of mean genetic diversity values between Algeciras and Tangier populations of Erica australis and Erica 
arborea. 







Erica australis 






Erica arborea 






Algeciras (N = 12) 


Tangier (N = 11) 


P-value 


Algeciras (N = 11) 


Tangier (N = 11) 


P-value 


A* 1 


8.683 


8.119 


0.103 


5.419 


5.212 


0.186 


Ho 2 


0.638 


0.618 


0.291 


0.576 


0.555 


0.213 


Hs 2 


0.754 


0.698 


O.OOI 


0.605 


0.589 


0.082 


F,s 


0.155 


0.116 


0.024 


0.048 


0.056 


0.758 


Fst 


0.057 


0.070 


0.278 


0.026 


0.038 


0.186 



M* = allelic richness calculated after the rarefaction method of El Mousadik and Petit (1996). 

2 H 0t H s , average values of observed heterozygosity and genetic diversity within populations, respectively. A/ - Number of sampled populations. Significant values based 
on 10,000 permutations are indicated in bold. 
doi:1 0.1 371 /journal.pone.0098602.t002 
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Figure 2. Bayesian analyses of genetic structure from twenty-three populations of £ australis (a, b) and twenty-two of £ arborea (c), 
respectively. Mean proportion membership of each population to the predefined, K=2 and K=4 (only E. australis) clusters with the highest AK 
values obtained following Evanno ef al. (2005). Proportion of membership of individuals for predefined clusters in both species. 
doi:1 0.1 371 /journal.pone.0098602.g002 
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Figure 3. Principal Coordinates Analyses (PCoA). It has been constructed with D A genetic distance of Nei et al. (1983), showing the association 
among populations of E. australis (a) and £ arborea (b). Black dots and white dots represent Algeciras and Tangier sites, respectively. Minimum 
Spanning Trees (MST) superimposed onto PCoA plots show the minimum length links between populations. 
doi:1 0.1 371 /journal.pone.0098602.g003 



subcluster was then linked to another, smaller group composed 
of the five southernmost Tangier populations (Fig. 3b). 

Populations of E. australis showed significant isolation by 
distance (IBD) within each side of the Strait of Gibraltar 
(Fig. 4a). Nonetheless, a stronger correlation between genetic 
distance (pairwise linearized F ST ) and geographical distance was 
obtained for Algeciras populations, as denoted by their signifi- 
candy steeper correlation line (/(-value <0.03, ANCOVA; Fig. 4a). 
In E. arborea, by contrast, a stronger correlation between genetic 
and geographical distances was obtained in Tangier (/rvalue < 
0.01, ANCOVA; Fig. 4b), whereas no significant correlation was 
obtained for Algeciras populations (Fig. 4b). 

Regarding estimated levels of gene flow between populations 
(jV m ), both E. australis and E. arborea showed conspicuous 
differences between the two sides of the Strait. In E. australis, a 
significant decrease in gene flow levels was found in Tangier 
between pairs of populations separated less than 30 km from each 
other (mean jV m , 95% CI: 5.08, 3.12-7.07) compared to Algeciras 
(9.73, 7.55-12.0), as evidenced by the non-overlapping confidence 
intervals. This difference vanished at distances longer than 30 km 
(3.86, 3.04-4.66 vs. 3.57, 2.96-4.11 in Tangier and Algeciras, 
respectively). In contrast, E. arborea showed similar levels of gene 
flow between populations at each side of the Strait, both for 
population pairs at distances shorter than 30 km (31.30, — 1.89— 



68.01 vs. 25.47, - 1.97-51.02) and longer than 30 km (7.92, 5.05- 
10.78 vs. 1 3.58, 6.96-20.24, in Tangier and Algeciras, respectively). 

Discussion 

Habitat patchiness leads to small population size and spatial 
isolation, which may also cause loss of genetic diversity in those 
populations by genetic drift and reduced gene flow among them 
[3,7]. This study illustrates the effect of an overall increase in 
habitat patchiness on the population genetic structure of habitat 
specialist species [17,56] compared to generalists [19]. It might 
seem contradictory, though, that the sandstone specialist E. 
australis has higher genetic diversity values than the generalist E. 
arborea (Table 2). This is likely a consequence of their different 
biogeographical histories and times of origin. The two species 
seem to be phylogenetically close [27,28], but E. arborea occupies a 
more derived position compared to E. australis in the published 
phytogenies [27,28] that, despite the absence of lineage dating 
analysis, suggests an earlier origin for E. australis. This species is 
endemic to the western end of the Mediterranean Basin and most 
likely originated somewhere in this region, on either side of the 
Strait of Gibraltar. In contrast, E. arborea is a widespread, circum- 
Mediterranean species [30], which probably originated in eastern 
Africa/ Arabia and reached the Strait area through recent 
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Figure 4. Isolation by distance (IBD). Matrix correlations between 
pairwise linearized F ST values (Slatkin 1995) (y-axis) and log-transformed 
pairwise geographical (x-axis) distance values for £ australis (a) and £ 
arborea. (b). Correlation values for £ australis within each side of were 
r = 0.654, p = 0.001 and r = 0.282, p = 0.046 for Algeciras and Tangier, 
respectively. Correlation values for £ arborea within each side of were 
r = 0.17, p = 0.142 and r = 0.475, p = 0.004 for Algeciras and Tangier 
respectively. P values reported after 1000 random permutations Mantel 
tests. Black circles and white circles indicate Algeciras and Tangier 
pairwise comparisons, respectively. Solid lines and dashed lines 
represent correlation lines for Algeciras and Tangier, respectively. 
doi:1 0.1 371 /journal.pone.0098602.g004 

(Pleistocene) range expansion [57]. The presumably later diver- 
gence of E. arborea and, particularly, a progressive loss of genetic 
diversity associated with its westward migration [57,58] would 
account for its overall low genetic diversity values in the Strait of 
Gibraltar compared to those of E. australis, both within and among 
populations (Tables 1 and 2). Nonetheless, what we highlight in 
this study is that average levels of within-population genetic 
diversity (H s ) and between-population gene flow levels (jV m ) in the 
sandstone specialist E. australis were significantly lower in Tangier 
(overall high sandstone patchiness) than in Algeciras (overall low 
patchiness), whereas no differences concerning sandstone patch- 
iness were detected in the edaphic generalist E. arborea (Table 2) 
between both sides of the Strait. 

The lower inbreeding coefficients (-PIS) in Tangier populations 
of E. australis than in Algeciras ones (Table 2) may at first seem 
contradictory, since inbreeding is expected to increase with 
fragmentation in plant populations [3]. However, _F IS may be 
biased upward in large, or highly spatially structured populations 
because of the Wahlund effect [3,59]. Therefore, F m should be 
considered with caution as an isolate indicator of inbreeding [7] 
and will not be discussed further. 

Low genetic diversity values in discontinuous populations are 
normally coupled with high levels of genetic differentiation [60]. 



This study did not detect significandy higher genetic differentia- 
tion between populations (F ST ) of the sandstone specialist E. 
australis in Tangier despite their significantly lower genetic diversity 
(H s ) values in this side of the Strait (Table 2). This may be partially 
a consequence of the mixed origin of populations in both sides of 
the Strait (Fig. 3a). Nonetheless, this species showed a somewhat 
higher genetic structure in Tangier (three genetic clusters are 
predominant in Tangier vs. two in Algeciras, Fig. 2b). By contrast, 
an increase in genetic structure in Tangier was not found for E. 
arborea (Fig. 2c), probably as a consequence of its broader 
ecological (edaphic) niche. Populations growing outside the 
sandstone patches would make population connectivity in this 
species more similar between both Tangier and Algeciras ranges 
than in the sandstone specialist E. australis. 

While significant correlations between geographical and genetic 
distances (i.e. significant IBD signals) were detected for the 
sandstone specialist E. australis at each side of the Strait (Fig. 4a), 
the IBD signal was only significant in Tangier for the generalist E. 
arborea (Fig. 4b). The wider edaphic spectrum of the latter species is 
one likely reason that might contribute to homogenize its genetic 
spatial structure under conditions of low edaphic patchiness. On 
the other hand, the presence of two genetic groups of populations 
of E. arborea in Tangier with a tendency of a north-south 
geographic structure (Fig. 2c; Fig. 3b) might account for the 
increase in IBD at this side of the Strait. Unlike E. arborea, the IBD 
signal somehow weakened in Tangier for E. australis, as reflected 
by the more gentle correlation slope comparisons of their 
correlation coefficients (Fig. 4a). This weakening of the IBD signal 
in Tangier populations of the sandstone specialist species is caused 
by strong genetic differentiation (F ST ) values between pairs of 
geographically close populations that fall in different, isolated 
sandstone patches (e.g. T12, T15, T16, T19, T23; see Fig. 1), 
indicating low levels of gene flow between the intervening 
populations [61]. Restrictions on gene flow are certainly a primary 
consequence of spatial isolation [14,62]. In this regard, the effect of 
the edaphic (sandstone) isolation in Tangier is clearly evidenced by 
the significant drop of relative gene flow levels (A" m ) between 
populations of the sandstone specialist E. australis distant less than 
30 km in this region compared to Algeciras, which was not found 
for the edaphic generalist E. arborea. 

High patchiness leads to genetic impoverishment and isolation 
of populations, which makes them more prone to (local) extinction 
[3,63]. These results provide evidence of how natural habitat 
(sandstone) patchiness has a more marked effect on the genetic 
diversity and gene flow in populations of a sandstone specialist 
heath species (E. australis) than in those of an edaphic generalist (E. 
arborea). Since most endemic species in Mediterranean heathlands 
of the Strait of Gibraltar are sandstone specialist [21], this study 
allows us to infer a high vulnerability of this endemic flora to 
genetic erosion and local extinction of populations in the Tangier 
side of the Strait of Gibraltar. 

Supporting Information 

Figure SI (a, b) Log-likelihood of the seven microsatellite loci 
data for 23 populations of E. australis (a) and of the eight 
microsatellite loci data for 22 populations of E. arborea (b) given K 
clusters, obtained through 10 runs of the STRUCTURE analysis 
for each K. Corresponding AK estimation (c, d) according to 
Evanno et al. (2005) showing maximum peaks of A^values at K= 2 
and K=\ for E. australis (c), and at K=2 for E. arborea (d), 
indicating that those are the optimal solutions for K given the data. 
(DOC) 
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Table SI Allele frequencies for 7 micro satellite loci in 

23 Erica australis populations. Numbers in brackets indicate 
sample sizes. Number of sampled individuals (JV) in each 
population is indicated in paretheses. 
(DOC) 

Table S2 Allele frequencies for 8 micro satellite loci in 
22 Erica arborea populations. Numbers in brackets indicate 
sample sizes. Number of sampled individuals (JV) in each 
population is indicated in paretheses. 
(DOC) 
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